Universal and measurable entanglement entropy in the spin-boson model 
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We study the entanglement between a qubit and its environment from the spin-boson model with 
Ohmic dissipation. Through a mapping to the anisotropic Kondo model, we derive the entropy of 
entanglement of the spin E(a, A,h), where a is the dissipation strength, A is the tunneling amplitude 
between qubit states, and h is the level asymmetry. For l-a> A/u c and (A, h) <C uo c , we show 
that the Kondo energy scale Tk controls the entanglement between the qubit and the bosonic 
environment (u c is a high-energy cutoff). For h <C Tk, the disentanglement proceeds as (/i/Tx) 2 ; 
for h 2> Tk, E vanishes as (Tk / h) 2 ~ 2a , up to a logarithmic correction. For a given h, the maximum 
entanglement occurs at a value of a which lies in the crossover regime h ~ Tk- We emphasize the 
possibility of measuring this entanglement using charge qubits subject to electromagnetic noise. 

PACS numbers: 03.65.Ud, 03.67.-a, 72.15.Qm, 85.3e.Be 



The concept of quantum entropy appears in multi- 
ple contexts, from black hole physics^ to quantum in- 
formation theory, where it measures the entanglement 
of quantum states^ Prompted by the link between en- 
tanglement and quantum criticality, 3 a number of re- 
searchers have begun to study the entanglement entropy 
of condensed matter systems. In this Letter, we em- 
ploy the spin-boson model 4 ^ to describe the entangle- 
ment between a qubit (two-level system) and an infinite 
collection of bosons. With an Ohmic bosonic bath, the 
spin-boson model undergoes a quantum phase transition 
of Kosterlitz-Thouless type when a — 1 = A/co c , where 
a is the strength of the coupling to the environment, 
A is the tunneling amplitude between the qubit states, 
and lo c A is an ultraviolet cutoff^! 7 - When the two 
levels of the qubit are degenerate, the entanglement be- 
tween the qubit and the bosons is discontinuous at this 
transition*^ Here we report the first rigorous analytical 
results for the entanglement (quantum entropy) in the 
strongly entangled regime l-a» A/w c . 

We exploit a mapping between the spin-boson model 
and the anisotropic Kondo model; our results follow from 
the Bethe ansatz solution of the equivalent interacting 
resonant level model. 10 ' 11 We show that the entropy of 
entanglement (E) of the qubit with the environment is 
controlled by the Kondo energy scale Tk, which gov- 
erns the low-energy regime of the Kondo problem (a 
strongly correlated Fermi liquid). We derive simple uni- 
versal scaling forms for the entanglement in the limits 
h -C T K and h > T K (Fig. 1), where h < lu c is the 
level asymmetry between qubit states. We also observe 
that, for a given h, E is maximized at a value of a 
which lies in the crossover regime h ~ Tk- While the 
spin-boson model describes many systems of experimen- 
tal interest^ the example most pertinent to this work is 
a noisy charge qubit, built out of Josephson junctions- 12 
or metallic islands ; n i 13 where the environment embodies 
the electromagnetic noise stemming from Ohmic resistors 
in the external circuit. 14 When the qubit and the leads 
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FIG. 1: (color online) Summary of our results. The shaded 
region depicts the crossover regime h ~ Tk- We compute 
Tk from the universal scaling function of Ref. I2H with A = 
0.1w c ; for the sake of clarity we choose a large value of A. 



form a ring, the entropy of entanglement can be con- 
structed from two measurable quantities: the persistent 
current in the rin g 11 ! 15 and the charge on the dotf^ 

Model and entanglement entropy. — The Hamiltonian 
for the spin-boson model with a level asymmetry h is: 
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where o~ x and a z are Pauli matrices and A is the tunnel- 
ing amplitude between the states with o~ z = ±1. H osc 
is the Hamiltonian of an infinite number of harmonic os- 
cillators with frequencies {u q }, which couple to the spin 
degree of freedom via the coupling constants {A 9 }. We 
assume an Ohmic heat bath with the spectral function 
J(uj) = irj^q^qtii^q — co) = 2iTauj, oj <C uj c . The di- 
mensionless parameter a measures the strength of the 
dissipation. For h = and A/u c <C 1, this model has a 
quantum critical line along the separatrix a — l = A/lu c & 
The region a — 1 > A/u c is a broken-symmetry phase 
(the "localized" phase) where A renormalizes to zero and 
lim; l ^o( cr z) 7^ 0; here the bosons disentangle from the 
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spin4 The "delocalized" phase (a — 1 < A/ui c ) is divided 
into two regimes by the separatrix 1 — a = A/uj c . The 
localized phase can be treated by perturbation theory in 
A, but in the delocalized phase this works only when h 
is larger 1 - We focus on the regime 1 — a > A/oj c , where 
the entanglement between the qubit and the environment 
leads to a renormalizcd tunneling amplitude A ron < A. 6 
At zero temperature, the entanglement between two 
members (A and B) of a bipartite system in the 
pure state \ip) is given by the von Neumann entropy 
E = -Trp A log 2 PA = -Trp B log 2 p B , where p A (B) = 
Tr B f A \\ip){ip\£l{ \ip) is the ground state of Hsb and A is 
the qubit, this results in E = —p + \og 2 p+ — p_log 2 p_, 

where p± = ^1 ± yj (a x ) 2 + (cr z ) 2 ^j /2; (a y ) = because 



Hsb is invariant under o~ y 



-Oy. We present exact 



results for E(a, A, h) in the regime 1 — a ^> A/tv c . Al- 
though E is defined at zero temperature, it exhibits uni- 
versality that is reminiscent of thermodynamic quantities 
like susceptibility and specific heat^ Recent work on the 
impurity entanglement in the isotropic Kondo model has 
emphasized universality in a similar vein.— 

Mapping onto the anisotropic Kondo model. — Our re- 
sults follow from a well-known mapping between Hsb 
and the anisotropic Kondo model^ defined as 

Harm = H kin + ^- ^ (c\^c k >iS~ + c) k[ c k ^S + 
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( c feT c fc'T ~ c ki c k'l) S z + hS z . (2) 
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This Hamiltonian describes the anisotropic exchange in- 
teraction between conduction electrons (labeled by the 
one-dimensional wave number k and spin a =],[) and 
a spin- 1/2 impurity. i?km is the kinetic energy of 
the electrons. Harm and Hsb are equivalent if we 
take A/u> c — * pJ±, a — ► (1 + 25/n) 2 , and h — ► h, 
where p is the density of states per spin of the elec- 
trons and 8 = t&n~ 1 (—irpJ z /4 : ) is the phase shift they 
acquire from scattering off the impurity.^ . The region 
1 — a > A/lo c corresponds to the antiferromagnetic 
Kondo model, while a — 1 > A/lu c corresponds to 
the ferromagnetic Kondo model. The equivalence be- 
tween Hs b and Hak m can be established via bosoniza- 
tion; a "re-fermionization" of Hsb then leads to an in- 
teracting resonant level Hamiltonian 20 which has been 
solved by Bethe ansatzJ^ The low-energy physics of the 
regime 1 — a 3> A/uj c is controlled by the Kondo scale 
T K = A(A/D) a ^ 1 - a '> ~ A rcn , where D is a high-energy 
cutoff. D and to c are different but their relationship is 
fixed — see, e.g., Ref. Oj]. Note that a more general ex- 
pression for Tk can be obtained from the renormaliza- 
tion group equations for a and A/w c j&2i in the limit 
1 — a — > A/oj c , Tk assumes the exponential form of the 
isotropic, antiferromagnetic Kondo model (Fig.[T]). 

Generalities. — It is clear from Eq. ([TJ) that (a x ) = 
—2dE g /dA and (a z ) = 2dE g /dh, where E g is the ground 
state energy of Hsb- Since Hsb and Hakm are related 



by a unitary transformation, they have the same ground 
state energy (up to an unimportant constant). The field 
h couples directly to the spin in both models, so we have 
(S z ) — (a z )/2. However, a similar relationship does not 
hold between (S x ) and (a x ), and therefore E does not 
measure the entanglement between the Kondo impurity 
and the conduction band. At a — 0, the qubit is decou- 
pled from the environment; thus, (a x ) — A/^fh? 
-h/Vh? 



+ A 2 
0, we 



and (a z ) = -h/\/Ii 2 + A 2 . With p+ = 1 and p- = 
have E = for all values of A and h. But when a 
for h — and A/ui c — + 0, the system is equivalent to 
the antiferromagnetic SU(2) Kondo model with J± w J z 
and (a x ) w (a z ) = 0, so we expect E — > 1, in agreement 
with previous Numerical Renormalization Group (NRG) 
results^ On the other hand, we must have E — > at 
large h, because the qubit is localized in the state with 
(a z ) = — 1 and (a x ) = 0. We argue that the Kondo map- 
ping allows us to examine how E interpolates between 
these limits and to explore the phase diagram (a, h). 

Toulouse limit. — First we focus on the point a = 1/2, 
which corresponds to the Toulouse limit of the Kondo 
model.— >22, The resonant level is non-interacting in this 
limits , so the ground state energy is simply that of a 
level at energy h with width ~ Tk- We find 
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First, consider the limit h <C Tk, where 
-{2/n){h/T K ) and (<r x 
\n(T K /D)}. The result for (cr 2 



-(4/7r)VWD[l + 
is consistent with the 
Kondo ground state, where S is fully screened and (S z ) cx 
h/TK at small h. Since both (o~ x ) and (a z ) are small, the 
system is close to maximal entanglement: 



lim £(1/2, A. h) 



E ^ 2 ^-^\t-k 



(5) 

where £(1/2, A, 0) = 1 - [l + ln(jk)] 2 . We 

have argued that E(h = 0) — > 1 as a — > 1 ; since the 
correction is already small at a = 1/2, we anticipate that 
E varies smoothly from a = 1/2 to a = 1 at h = 0. The 
second term in Eq. ([5]) is a universal function of H/Tk, 
with a quadratic dependence on energy that arises from 
the Kondo Fermi liquid behavior of (a z ). In the opposite 
limit h > T K , we find that (a z ) -> -1 + 2T K /{ith) and 
(a x ) -> -(4/Tv)y/T K /Dln(h/D). Again, the leading h- 
dependence of E has a universal form dictated by (a z ): 
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Because of the logarithmic correction, E approaches zero 
slowly at large h. Plots of E over the full range of h are 
shown in Fig. [2] for several values of Tk- 
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FIG. 3: The coefficients fci(a) (left) and k 2 (a) (right). 



FIG. 2: (color online) Entropy E(a = 1/2, h), plotted on a 
logarithmic scale for five values of Tk = A 2 / D; from top 
to bottom T K /D = 0.005, 0.001, 0.0005, 0.0001, 0.00005. The 
solid lines show the asymptotes found in Eqs. ((5| and (|6} . 



Away from a — 1/2. — The Bethe ansatz solution 
of the interacting resonant level model provides exact 
solutions for (a z ) and (a x ) in the delocalized realm 
1 — a 3> A/icic ! 10 ! 11 While the general expressions are 
quite complicated, we derive simple scaling forms for the 
entanglement entropy in the limits h <C Tk and h 3> Tk- 

For /i « Tk wc find 



lim (a z 

h<£T K 



2 e 2(T^o p[l + 1/(2 - 2a)] / h 



r[l + a/(2 - 2a)] 



where b = a In a + (1 — a)ln(l — a). Again we have 
(a z ) oc h/Tx at small h, in keeping with the Kondo 
ground state. The leading /i-dependence of (a x ) is of 
order h 2 and therefore negligible, which leaves 



lim (a x ) = 1 — + Ci(a)-^-, 

ft^CTjf 2a - 1 w c A 
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with C\ (a) 
T K 



_ e -b/(2-2a) r[l-l/(2-2a)] 



As 



0F(l-a) r[l-a/(2-2a)] ■ u ^ 0, 

A and Ci(0) = 1, so we recover the exact result 
(c x ) a= h=o — 1 U P to a correction of order A/uj c . This 
ensures E — > for a — > 0. As we turn on the coupling to 
the environment, we introduce some uncertainty in the 
direction of the spin and (a x ) progressively decreases. 
For a < 1/2, the monotonic decrease of Tk/A domi- 
nates. For a > 1/2, the first term in Eq. ([5]) dominates 
and we have (a x ) ~ A/lo c w 0. The smallness of (a x ) in 
this regime reflects the loss of coherent Rabi oscillations 5 
that occurs at the dynamical crossover a = 1/2. Note 
that (a x ) remains analytic: in the limit a — > 1/2, we 
take Ci(a) = (4/vr)r(l - 2a) 4/(vr(l - 2a)) and 
use the identity D(a = 1/2) = 4oj c /tt to find (<r x ) — > 
- ( 4 A) V t k /Dln(T K /D), in agreement with Eq. ©. 

Now we focus on the region away from a = 0, where 
Tk <C A and the system is strongly entangled at h = 0. 
Here we can generalize Eq. ([5]): 
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The coefficient k\ (a) (Fig. [3]) is given by 
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and E(a, A,0) = ^ - 21E2 V^T^I 
{Ji/Tk) 2 scaling, which is a feature of the Fermi liquid 
fixed point, persists for a ^ 1/2. Note that the scaling 
with h/TK is determined entirely by (cr z ), while the (non- 
universal) contribution at h — arises from (o~ x ). We also 
observe that E(h — 0) saturates at maximum entropy 
for a > 1/2, where the leading correction is of order 
(A/oj c ) 2 . The plateau for a > 1/2 demonstrates a link 
between entanglement entropy and decoherence. In the 
limit a — > 0, the leading ^.-dependence is still quadratic, 
but with a non-universal pre-factor. 
For h 3> Tk we have 
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where 62(0) 



r(3/2-a) 



In the limit a — ► 1/2, 



ViF(l-2ce) r(l-a) 

(ctj;) contains two additional terms which conspire with 
the other terms to produce the logarithm of Eq. (4); we 
do not write them explicitly because they cancel each 
other for a 7^ 1/2. In the regime Tk <C h <C A we find 



lim E(a. A, h) = fc 2 (a) In , 
T if «/ t «A v ' ' ; 2y ' \T K 



where the pre-factor (Fig. [3J is given by 
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As in Eq. ([6]) , the universal scaling function follows from 
the high-field response of (<J Z ). 

In the localized phase we obtain (a z ) « —1 and 
(a x ) w — and therefore -E w — for infinitesimal /i. 
Since dissipation localizes the spin in the " down" state for 
h = + , we do not expect the entropy to depend strongly 
on external field; so it makes sense that k 2 — ► as a — > 1. 
Deep in the localized phase, for a — 1 ^> A/uj c , we find 
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FIG. 4: (color online) E(a, A = 0.01a;c, h) versus a at several 
values of h. At h = 0, we check that i? cx a in the limit a — > 0. 
The arrow marks the value of a at which Tk = 0.001lj c ; we 
see that for h = 0.001cj c , E is maximized near this point. 

from perturbation theory-^ that E ~ — (A/o; c ) 2 ln(A/w c ) 
to leading order for all h (Fig. 1). As we approach the 
critical line a — 1 = A/lu c from the localized side, this 
behavior is replaced by E ~ — (A/w c ) ln(A/w c )i^ 

We use the full solution o 10 ' 11 for (err) and (cr z ) to plot 
.E versus a for various /i in Fig. [5] The entropy increases 
monotonically when h = 0: it is linear in a near a = 0, 
and it saturates at £ ~ 1 for a > 1/2, as discussed 
above. As h increases from zero, E exhibits a maximum 
at progressively smaller values of cv, in agreement with 
previous NRG results.— In our view, this maximum signi- 
fies the crossover h ~ Tk (Fig- El- If h — 0, the entropy 
E is driven to zero by dissipation, and we observe instead 
a sharp non-analyticity at the phase transition. 8 

Experiments. — An important open question in the 
study of quantum entanglement is whether it can be mea- 
sured experimentally. The model considered here is re- 
alized in noisy charge qubits, composed either of a large 
metallic do t 11 ! 13 (the single electron box) or a supercon- 
ducting island^ (the Cooper pair box). The gate volt- 
age controls the level asymmetry h, and A corresponds 



to the tunneling amplitude between the dot and the lead 
or the Josephson coupling energy of the junction. If the 
gate voltage source is placed in series with an external 
impedance, voltage fluctuations will give rise to dissipa- 
tion even at zero temperature^ The parameter a can 
be varied in situ when a two-dimensional electron gas 
acts as the Ohmic dissipative environment^ Here, E 
depends only on (a x ) and (cr z ), so it can be constructed 
from physical observables. While these quantities would 
obviously be measured at finite temperature, we assume 
it is possible to recover the ground state density matrix 
by extrapolating them to their zero-temperature values. 
Charge measurements^ yield the quantity (o~ z ), which 
represents the occupation of the dot or island. In a ring 
geometry, the application of a magnetic flux generates a 
persistent current that is proportional to the observable 
(a x )£L Another promising system is the atomic quan- 
tum dot, which also permits experimental control of the 
coupling between the dot and the bosonic reservoir^ 

Conclusion. — We have provided quantitative predic- 
tions for the entropy of entanglement of the spin in the 
spin-boson model. This entropy exhibits universal behav- 
ior in the delocalized phase, governed by the Fermi liquid 
fixed point of the equivalent anisotropic Kondo system. 
We have also described an experimental setup capable of 
testing our predictions; such measurements would pro- 
vide an empirical proof of the existence of entanglement 
entropy. Although the presence of dissipation in charge 
qubits makes them unlikely candidates for a functioning 
quantum computer, we have shown that they can be used 
to explore interesting links between quantum entangle- 
ment, decoherence, and quantum phase transitions. This 
work might be extended to two noisy qubits. 
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